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^ | Abstract 

^ ■ We discuss questions related to the application of OPE for the long-lived systems with 

q . two heavy flavors. The values of quark masses entering such calulations are constrained by 

use of currently available data on the lifetimes of and 5 c -mesons. The lifetimes of doubly 
■ charmed baryons at these values of parameters are evaluated. The necessary comments on the 

difference in the results obtained previously are made. 

> : 

ctj ! 1 Introduction 

The main motivation for studying the weak decays of doubly heavy systems is to test our under- 
CT\ standing of QCD in the limit, where some masses involved are heavy and, as a consequence, certain 
aspects of bound state dynamics are simplified. In the case, when the decay of bound system pro- 
ceeds due to electro-weak interactions, the consideration also gives us a possibility to extract some 
basic properties of quark interactions at a fundamental level, including precise determination of 
CKM parameters. The analysis of decays for the doubly heavy systems already has a long history. 
^ First of all, the decays of J I ip and T were considered. Being composed by the quark and antiquark 
>• . of the same flavor, these mesons decay mainly trough quark- antiquark annihilation into hadrons or 
lepton pair. Another particle, acquiring a lot of attention, was the 5 c -meson It represents a first 
long-lived particle in the family of doubly heavy systems, whose decays proceed due to the weak 
interactions. The other representatives in the family of doubly heavy hadrons are baryons with two 
heavy quark, yet to be discovered experimentally. In this paper we would like to consider lifetimes 
and some issues of OPE, used in such the framework for the hadrons with two heavy quarks, which 
decay due to the weak interactions. 

Weak decays of the ground state of £? c -meson together with semileptonic and various exclusive 
modes were considered in |3|, || [7], ^, [!|. The estimates of total £? c -lifetime in various quark 
models were done in j|, |^, |6| . The first OPE based result for the lifetime of B c - meson was obtained 
in [|T(J , where, however, the matrix element for the dimension 5 operator Qga ■ GQ was incorrectly 
evaluated, which led to the overestimation of 5 c -lifetimeQ. In addition, the leading corrections 
to the spectator decays of b and c-quarks determined by the kinetic energy, weak annihilation and 
Pauli interference, were not calculated explicitly in [|K|. A systematical OPE-based approach to the 
evaluation of 5 c -meson lifetime was developed in ||11||. In a short time, the lifetime value predicted in 



1 See a discussion of this point in [IT] 
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TTJ was experimentally confirmed by the discovery of 5 c -meson at FNAL by the CDF-collaboration 



EJ- 

At the next stage in studying the decays of doubly heavy systems the approach of fll 
generalized to the case of baryons with two heavy quarks in 
at different values of parameters was done in |I4[ 



was 
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A repetition of our results 
The questions related to the spectroscopy, sum 

corresp ondingly. 
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rules and production of doubly heavy baryons were considered in 

In this work we discuss the question related to the choice of parameters used in the OPE-based 
frameworks for the doubly heavy systems, as required to obtain the theoretical predictions on the 
lifetimes in agreement with the experimentally measured values. 

This paper is organized as follows. In Section 2, we define notations in the framework of OPE- 
calculations for the total widths of hadrons with two heavy quarks. In Section 3, the procedure 
for estimatinf the non-perturbative matrix elements between the states with two heavy flavors is 
considered. Section 4 is devoted to the numerical evaluation and discussion of parameter choice 
explored for the calculation of lifetimes for the doubly charmed baryons to be consistent with 
the experimentally measured value of 5 c -meson lifetime. We draw conclusions in Section 5 by 
summarizing our results. 



2 The OPE framework for lifetimes 

In accordance with the optical theorem, the total width Th for the hadron H, where H is the 
B c - meson or one of the baryons and Q cc , has the form 

1 



H 



2M 



■(H\T\H), 
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where we accept the ordinary relativistic normalization of state, (H\H) = 2EV, and the transition 
operator T: 



T = 3m / d 4 x {TH eff (x)H eff (0)} 



(2) 



is determined by the effective lagrangian of weak interaction H e ff at the characteristic hadron 
energies: 

Heff = ^=V q2qs V* qi [CMO + + CL(/*)CL] + h.c. 

where 
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so that / denotes the number of flavors, and Q marks the flavor of heavy quark (b or c). 

The quantity T in ([!]) permits the Operator Product Expansion in the inverse powers of heavy 
quark mass. The reason is that, the energy release in the weak decay of either quark is large 
compared to the scale of bound state dynamics, and, so, we can expand in series over the ratio of 
these scales. In this way, the OPE has the form: 



i=l 



m 



l —CMQ l go^Qt + ^O(l)}. 
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The leading contribution is given by the spectator decay, i.e. by the term QQ, which is the 
operator of dimension 3. The corrections to the spectator decays of quarks are given by the operator 
of dimension 5: Qgq = Qg&nuG^Q and operators of dimension 6, Q2Q2q = QTqqT Q, where the 
dominant contributions are provided by the Pauli interference and weak scattering. For the latter, 
we have 

T B + = T 35b + T 35c + T 6 ^j + T$f A , 

where T 35 q denotes the contributions into the decays of quark Q by the operators with the di- 
mensions 3 and 5, and the forthcoming terms are the interference, scattering and annihilation of 
constituents. In the explicit form we find 

T 35b = V hsp Jb - %2P cl + P CTl + K ob (P cl + P ccl ) + K 2b (P c2 + P cc2 )]0 Gb} (5) 

% 5c = r c , spec cc - %(2 + K 0c )P sl + K 2c P s2 }0 Gc , (6) 
m c 

where 

with K q = C 2 _ + 2C 2 ., K 2Q = 2(C 2 _ — C 2 ), and T QtSpec denotes the spectator width (see |1| [2~lj 
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P cl = {l-y)\ P c2 = (l-yf, 



Peri = v / l-2(r + y) + (r-y) 2 [l-3(r + y) + 3(r 2 + y 2 )-r 6 -y d -4ry + 

7ry(r + y)] + 12rV In (1 -r -y + ^ -2(r + y) + (r -y)^ (g) 



Pea = v 7 !-^! - Qy + 2y 2 + 12y 3 )2% 4 In j-± (10) 



Pcc2 = \A - 4y(l + J + 3?/ 2 ) - 33/(1 - 2y 2 ) In = ' v (11) 
2 1 - VI - % 

2 

where y = ^§ and r = m 2 /m 2 . The functions P s i{P S 2) can be obtained from P c i(P C 2) by the 

substitution y = m 2 /m 2 . In the 6-quark decays, we neglect the value m 2 /m 2 and suppose m s = 0. 

The calculation of Pauli interference for the products of heavy quark decays with the quarks in 
the initial state, weak scattering and annihilation of quarks, composing the hadron, results in: 

^6,PI = ^PI,cs (12) 
^6,WA = T W a,cs + TwA,ud + %VA,vil (13) 
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T&i = ?Tpim ( 14 ) 

(15) 

%% = ^PiM+ 2 Y. r Pi^i ( 16 ) 



so that 



27^- = - *' v/ ' m(l--)K + m„)10(l-^--(l- Z _r x 

((C+ + C_f + ^(1 - 4£^)(5C 2 + C 2 _ - 6C_C+)) (17) 
^ 3G||^(0)| 2 2/ m dx2/ w x2 

(c 2 +ci + i(i-4^)(^-c 2 ; 

- 13G 2 F |^(0)| 2 2/1 m s , 2/ 

2T£ T s = - ^^m 2 (l V + m s x 

pi,ud Un cV m J V 

((C + - C_) 2 + ^(1 - 4A;5)(5C 2 + C 2 + 6C_C + )) (19) 

G\ |^(0)| 2 9/ m SN9/ , r . x9 17, ,0^1 . . 

m c I 3 J 

Tj 7iCS - = ^|K b | 2 /lM Bc K-m c ) 2 (l-^) 2 (2C 2 -C 2 _) (21) 



= 



12tt' 
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= ^|K b | 2 /l c M B m 2 (l-^ + ) 2 (4C 2 +C 2 +4C + C_) (22) 



= ^\V cb \ 2 f B M Bc m 2 T {l-z T ) 2 (23) 

Tpi, Ue e = ^pi^fi = ^~pi,u t t( z t — 0) (24) 

Here ^(0) is the value of quark-diquark baryon wavefunction at the origin, and /b c is the leptonic 
constant for 5 c -meson. In the evolution of coefficients C + and C_, we have taken into account the 
threshold effects, connected to the heavy quark masses. 

In expressions (|J) and @, the scale /x has been taken approximately equal to m c . In the Pauli 
interference term, we suggest that the scale can be determined on the basis of the agreement of 
the experimentally known difference between the lifetimes of A c , and S° with the theoretical 
predictions in the framework described aboveP]. In any case, the choice of the normalization scale 
leads to uncertainties in the final results. 

At present, the spectator decays of heavy quarks, contributing to T35Q, are known in the logarith- 
mic approximation of QCD to the second order [^, [27], pSfl , including the mass corrections 



in the final state with the charmed quark and r-lepton [fjl| in the decays of 6-quark and with the 



strange quark mass for the decays of c-quark. In the numerical estimates, we include these correc- 
tions and mass effects, but we neglect the decay modes suppressed by the Cabibbo angle, and also 
the strange quark mass effects in the 6-decays. 



2 A more expanded description is presented in p 
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Thus, the calculation of lifetimes for the doubly heavy systems under consideraton is reduced 
to the problem of evaluating the matrix elements of operators, which is the subject of next section. 



3 Hadronic matrix elements 

By use of motion equations, the matrix element of operator QiQi can be expanded in series over 
powers of 1/rriQj: 



WWIHU = 1 - {HmW? - t G)l<y|g> - + 0(4-). (25) 

2ml,- mz 



Thus, we have to estimate the matrix elements of operators from the following list only: 

Q j (iD) 2 Q j , {\)Q j aGQ\ Q> la (l - i 5 )Q'qi a ^ ~ 7s)?, 

Q J ',r:,Q'<r n i\ - 7s)g, Q 3 iai,Q 3 Q k i a {i - i,)Q\ (26) 

The meaning of each term in the above list was already discussed in the previous papers on the 



decays of doubly heavy baryons |12], [L3[ and in |fLl]| , so we omit it here. 

Further, employing the NRQCD expansion of operators QQ and Qgcr^G^Q, we have 

QQ = *Q*Q-^h( iD T*Q + ^U iD ) 4 *Q- 

^ Q g*B* Q -^ Q (DgE)* Q + ... (27) 

Qga^GTQ = -2tf Q gvB* Q - — & Q (DgE)* Q + ... (28) 

Here the factorization at scale \i (mg > \i > ttiqVq) is supposed. We have omitted the term of 
tyq(r{gE x D)^q, corresponding to the spin-orbital interactions, which are not essential for the basic 
state of hadrons under consideration. The field \I/q has standard non-relativistic normalization. 

Further, the phenomenological experience in the potential quark models shows, that the kinetic 
energy of quarks practically does not depend on the quark contents of system, and it is determined by 
the color structure of state. So, we suppose that the kinetic energy is equal to T = m b v b /2 + m c v 2 /2 
in the 5 c -meson and to T = m^v J/2 + rriivf/2 in the case of doubly heavy baryons for the quark- 
diquark system, and it is T/2 — m b v 2 /2 + m c v 2 /2 in the diquark (the color factor of 1/2). Then 



(B c \%(iD) 2 V b \B c ) 2 2m c T 



2M B m\ b m b (m c + m b ) 



(29) 



(BM(iDf%\B c ) ^ v , ^ 2m b T (3Q) 



2Mb c v(i 2 c ' m c (m c + m b ) 



2M B * QQi m 2 Q (: {m q + m Q > + m Q )(m Q > + m Q ) m Q (m Q + m Q >)' 
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(E» QQI \¥ QI (iDy* Q ,\E* QQI ) 2m^T + m Q T ^ 



2M s ^ Q m 2 Q , Q ' (m q + mcy + m Q )(m Q i + m Q ) m Q ,{m Q + m Q >)' 

Applying the quark- diquark approximation for the doubly heavy baryons and relating the matrix 
element of chromomagnetic interaction of heavy quarks in the S c -meson and that of diquark with 
the light quark to the mass difference between the exited and ground states Mh* — Mh, we have 

{B c \cc\B c ) 1 , 3 M B * — Mb c / m b \ 

2M Bc 2 Vc 4 m c \ 2mJ 

« 1 - 0.190 + 0.037 - 0.061 + ... (33) 

(EllcclEl) 1 1 2 1M™. -M S o 5# 2 lTd/oXl2 
x ccl — ! — ££1 = 1 of =££. _ — \\SJ d (Q) 4- 

2M=. 2 c 3 m c 18m? 1 1 ;| ^ " 

'— 'cc L C 

« 1 - 0.073 - 0.025 - 0.009 + . . . (34) 
Mc« c > _ 1 l., 2 l Mngg-Mng, 5fl 2 d | ^ 



2M Q o 2 c 3 m c 18m3' 

« 1 - 0.078 - 0.025 -0.009 + ... (35) 



Our presentation here is less detailed than in previous papers [12], [L3|. However, we hope, that the 
interested reader can find there all needed details. Numerically, we have assigned T ~ 0.4 GeV. 
The values of |^ d (0)| and M H are given in the next section. 

Analogous expressions can be obtained for the matrix elements of operator Qgu^G^Q 

iE^pm = 3mc (M, - M B J (! - „ _ , 16 , (36) 

<2 



(S^cyG^cl^) = 4 (M sg * - MgoJ 7g 2 



^ d (0)| 2 « -0.124, (37) 



(O qq ,|c^G^c|O qq ,) = {~° QQI \cga, v G^c\~? QQI ) (38) 

The permutations of quark masses lead to the required expressions for the operators of bb and 
bga^G^b. 



4 Numerical results and Discussion 



The analysis of S c -meson lifetime and lifetimes of doubly heavy baryons considered in [II], [L2|, [13 
shows a strong dependence on quark mass values. As has been said in Introduction, the aim of this 
work is to reduce this uncertainty in a way for the case of doubly heavy baryons. For this purpose, 
we would like to note, that the OPE framework used by us for the calculation of lifetimes of doubly 
heavy baryons is a generalization of what was previously developed for theS c -meson lifetime ||11|| . 
So we can test different sets of parameters used in the calculations of doubly heavy baryons lifetimes 
to the case of 5 c -meson, for which we already have some experimental data. 

The total width of _B c -meson consists of the spectator b and c-quark decays, Pauli interference 
and weak annihilation contributions. The lifetime dependence on the 6-quark mass can be eliminated 
by the requirement, that for any value of m c , the m& mass value is obtained by the matching which 
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results in the B d - meson lifetime to be equal to the experimentally measured value r Bd ~ 1.55 ps. 
This prescription leads to the following approximate relation between the heavy quark masses 



m b = m c + 3.5GeV (39) 

The performed analysis shows, that the 5 c -meson lifetime dependence on the c-quark mass is quite 
large, and it is hard to eliminate the dependence by use of some external input. For example, a 
quite large c-quark mass or too low virtualityf] are required to reproduce the absolute lifetime and 
semileptonic width in the OPE-based approach. Note also, that the application of OPE is based on 
the assumption about the quark-hadron duality. The latter can be violated in the case of D-mesons 
2"9"| , |3U|] , and, thus, in this case the OPE predictions are less reliable. On the other hand, as was 



noted by authors of ||11|| , there is no obvious violation of this duality in the case of B c - mesons. 
So, here we attempt to extract the c-quark mass value by fitting the OPE result for the S c -meson 
lifetime to the experimentally measured value to use in the calculations of lifetimes for the doubly 
heavy baryons, yet to be discovered. 

In the calculations of 5 c -meson lifetime we have used the following set of parameters 

m s = 0.2GeV \V cb \ = 0.04 M Bc = 6.26GeV 

M B * - M Bc = 0.073GeV T = 0.37GeV f Bc = 0.5GeV (40) 

It is just the set of parameters, implemented in the original paper on the B c lifetime. Here, we 
would like to note, that the quoted value of f Bc is bigger, than one obtained in the framework of 



QCD sum rules ]3T|, [32], [33J. 

In the calculation of quark spectator decays we put the renormalization scale to /i = nig and 
in the case of nonspectator decays we have Q /i = ni red , being the reduced mass for the S c -meson 
system. In Table [I] we have collected the numerical values of 5 c -meson lifetimes and relative 
spectator and nonspectator contributions for the different sets of b and c-quark masses, satisfying 



relation (|39|). Analyzing this Table, we see that, for example, the set of parameters, proposed in [I 



Parameters, GeV 


J2 b -> c, ps 1 


£c -> s, ps 1 


PI, ps" 1 


WA, ps" 1 


tb c , ps 


nib = 5.0, m c = 1.5, m s = 0.20 


0.694 


1.148 


-0.115 


0.193 


0.54 


nib = 4.8, m c = 1.35, m s = 0.15 


0.576 


0.725 


-0.132 


0.168 


0.75 


nib — 5.1, m c = 1.6, m s = 0.45 


0.635 


1.033 


-0.101 


0.210 


0.55 


nib — 5.1, m c = 1.6, m s = 0.20 


0.626 


1.605 


-0.101 


0.210 


0.43 


nib — 5.05, m c = 1.55, m s = 0.20 


0.623 


1.323 


-0.107 


0.201 


0.48 


nib — 5.0, m c = 1.5, m s = 0.15 


0.620 


1.204 


-0.114 


0.193 


0.53 



Table 1: The value of B c - meson lifetime together with the spectator and nonspectator contributions 
to the width at various choices of parameters. 

for the calculation of lifetimes for the doubly charmed baryons, is completely inconsistent with the 
3 At low values, the perturbation theory cannot be justified. 

4 As in the case of £? c -meson, the nonspectator contributions are not so large as compared to those of doubly 
charmed baryons, we have not introduced an additional low energy logarifmic renormalization for them. 
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experimental data, when applied to the calculation of i? c -lifetime. At this set, we have tb c =0.75 
ps contrary to 



T 



, ,, = 0.46 ± 0.18 s * a * ± 0.03 ss/s *ps. (41) 

The best set of parameters turns out to be = 5.05, m c = 1.55, m s = 0.2 GeV. So, now we in the 
situation, when our set of parameters has a strong motivation, and we may, at a confidence level, 
apply it to the calculation of lifetimes for the doubly heavy baryons. 
Estimating the lifetimes of doubly charmed baryons we put 

m s = 0.2 GeV m* = 0.45 GeV m* = 0.3 GeV | V cs | = 0.9745 

M =++ = M =+ = 3.478 GeV M Qcc = 3.578 GeV (42) 
Mbj* - M H?c = M n * c - M Qcc = 0.132 GeV T = 0.4 GeV ^ d (0) = 0.150 GeVi 

The numerical values of parameters, characterizing the doubly charmed baryons, are taken from 



IB] , |l~6|j . m* and m\ are the strange and light quark constituent masses, used for the calculation of 
bound state effects in the hadronic matrix elements of doubly charmed baryons. For the value of 
light quark- diquark function we assume 

1^.(0)1' = (43) 

where fjj — 170 MeV. This expression was obtained by performing the steps similar to |34], |3BJ for 
the derivation of hyperfine splitting in the light quark-diquark system. The factor k~9 accounts 
for the low energy logarithmic renormalization of fo constant. In the calculation of nonspectator 
effects, we have accounted for the low energy logarithmic renormalization to a hadronic scale //, 
where the latter was determined from the fit of theoretical predictions for the lifetime differences of 
baryons A c , and H° to the corresponding experimental values. Tables @,[|f| contain the values 
of lifetimes and relative spectator and nonspectator contributions for the doubly charmed baryons, 
as calculated at different sets of parameters used previously |I2], |i~4|| , and that of obtained from the 
fit of 5 c -meson lifetime to the experimentally measured value. 



Parameters, GeV 


Y^c^ s, ps" 1 


PI, ps" 1 


T = ++, ps 

—rr. 


m c = 1.35, m s = 0.15 


1.638 


-0.616 


0.99 


m c = 1.6, m s = 0.45 


2.397 


-0.560 


0.56 


m c = 1.55, m s = 0.2 


3.104 


-0.874 


0.45 



Table 2: The value of S+ + lifetime together with the spectator and nonspectator contributions at 
various values of parameters. 

The lifetimes, calculated at m c = 1.6 GeV, m s = 0.45 GeV, differ from those calculated previ- 
ously at these values of parameters, because of different value for the wavefunction of light quark- 



diquark system, used previously. The lifetime of differs from that of calculated in []rjj because 
of the weak scattering contribution to the total lifetime of baryon under consideration, which was 
wrongly estimated before in ||14j| . The final comment concerns with the importance of Pauli inter- 
ference in semileptonic inclusive decays of c-quark, which was introduced by Voloshin. As argued 
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Parameters, GeV 


s, ps" 1 


WS, ps" 1 


T~+ , ps 


m c = 1.35, m s = 0.15 


1.638 


1.297 


0.34 


m c = 1.6, m s = 0.45 


2.397 


2.563 


0.20 


m c = 1.55, m s = 0.2 


3.104 


1.776 


0.20 



Table 3: The value of lifetime together with the spectator and nonspectator contributions at 
various values of parameters. 



Parameters, GeV 


Y,c -> s, ps" 1 


PI, ps" 1 




m c = 1.35, m s = 0.15 


1.638 


1.780 


0.30 


m c = 1.6, m s = 0.45 


2.397 


0.506 


0.34 


m c = 1.55, m s = 0.2 


3.104 


1.077 


0.24 



Table 4: The value of f2+ lifetime together with the spectator and nonspectator contributions at 
various values of parameters. 



by the authors of [14], this term can be valuable. So, it is given by the following equation: 



T Voioshin = ^L|V r cd | 2 m^(4V^— l)5|^ ca (0)| 2 . (44) 

127T 

This term being doubly Cabbibo suppressed in the case of S+-baryons, does not give any sizeable 
contribution, when we discuss the lifetimes of these baryons and should be taken into account only 
in estimations of semileptonic branching ratios of heavy baryons. In the case of f2 cc -baryon it is no 
longer suppressed, and, so, this term is explicitly accounted for in our formulae. 

As can be seen from the previously performed analysis fll2 , p~3| , |14|| , the lifetimes of doubly 



heavy baryons strongly depend on the value of wavefunction for the light quark-diquark system at 
the origin. Also, we can present some arguments for its determination, but in all the cases these 
arguments rest on the hypothesis of quark-diquark picture for the doubly heavy baryons. In Fig. [I], 
0, |3| we have plotted the dependence of lifetimes for the doubly charmed baryons on this parameter. 
The precise value of light quark-diquark wavefunction is under question, and, thus, the uncertainty 
in its value should be included in the presented estimates of lifetimes. 

Finally we would like to comment on the theoretical errors in the given results. They are mainly 
caused by the following: 

1) The uncertainty in c-quark mass, taking into account the experimental errors on the B c - meson 
lifetime, can lead to jr ~ 15%. 

2) The uncertainty in the value of light quark-diquark wave function at the origin can lead to 
f«30%. 

Combining these two sources of uncertainties we get the total error of presented estimates at 
the level of 45%. 
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Figure 1: The dependence of -baryon lifetime on the value of wavefunction of light quark-diquark 
system at the origin |\l/ d '(0)|. 
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Figure 2: The dependence of S+-baryon lifetime on the value of wavefunction of light quark-diquark 
system at the origin \^ dl ' 
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Figure 3: The dependence of fi cc -baryon lifetime on the value of wavefunction of light quark-diquark 
system at the origin \ty dl < 



5 Conclusion 

In the present paper we have performed a detail investigation of parameter influence on the lifetimes 
of systems with two heavy flavors, calculated in the OPE-based approach. The fit to currently 
available data on the lifetimes of and 5 c -mesons, allowed us significantly to constraint the region 
of heavy quark masses. We present the numerical estimates for the lifetimes of doubly charmed 
baryons and discuss the huge difference between the results obtained previously. 
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